Goos-Hanchen shift and localization of optical modes in deformed microcavities 
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Recently, an interesting phenomenon of spatial localization of optical modes along periodic ray trajectories 
near avoided resonance crossings has been observed [J. Wiersig, Phys. Rev. Lett. 97, 253901 (2006)]. For the 
case of a microdisk cavity with elliptical cross section we use the Husimi function to analyse this localization 
in phase space. Moreover, we present a semiclassical explanation of this phenomenon in terms of the Goos- 
Hanchen shift which works very well even deep in the wave regime. This semiclassical correction to the ray 
dynamics modifies the phase space structure such that modes can localize either on stable islands or along 
unstable periodic ray trajectories. 
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I. INTRODUCTION 

The phenomenon of avoided level crossings has been dis- 
covered more than 70 years ago by J. von Neumann and E. P. 
Wigner IH]. Two energy levels undergo an avoided cross- 
ing under variation of a control parameter when they first ap- 
proach very close to one another, and then depart from each 
other without crossing. The energies Ej of eigenstates in 
closed quantum-mechanical systems are real-valued. In open 
systems, the energy eigenstates are replaced by quasi-bound 
or resonant states with complex energies Ej 121 [sl]. The imag- 
inary part determines the lifetime r, oc 1 /\m{Ej ) of the de- 
caying state. The corresponding generalization of avoided 
level crossings are avoided resonance crossings (ARCs) il. 
ARCs have been studied in a number of physical settings, 
such as laser-assisted electron-atom scattering |5], predisso- 
ciation dynamics of molecules |6], biased multiple quantum 
wells [7J, coupled quantum dots |8], quantum-dot microcavity 
systems IH 



llLmicrowave cavities lfl2i [Tsll . and optical 



microcavities |14i 115(1 . In the latter systems the electromag- 
netic modes and their frequencies play the role of eigenstates 
and their energies. One interesting phenomenon that can oc- 
cur near ARCs is the formation of fast and slowly decaying 
states [12., 161. 

In a recent letter it has been shown that this formation of 
fast and slowly decaying states can be accompanied by a lo- 
calization of the states along unstable (or marginally unstable) 
periodic trajectories of the underlying classical system ifisll . 
This localization is somewhat similar to the well-known effect 
of scarring originally discovered for closed chaotic systems in 
the field of quantum chaos 117, 18]. 

The aim of the present paper is twofold. First, we discuss 
the change of the phase space structure of modes in optical 
microcavities when such localization near an ARC occurs. 
The phase-space picture allows to clearly identify the physical 
mechanism of lifetime enhancement (reduction): destructive 
(constructive) interference in real space leads to a reduction 
(enhancement) of intensity in the leaky region of phase space. 
Second, we present a semiclassical explanation of the appear- 
ance of localized modes near ARCs based on an augmented 



ray dynamics in which the Goos-Hanchen shift (GHS) is in- 
cluded. The GHS is a lateral shift of totally reflected beams 
along the optical interface 119], i.e. the points of incidence 
and reflection do not coincide. Rather, the beam appears to 
travel, as an evanescent wave, a short distance through the 
optically thinner medium. The GHS is proportional to the 
wavelength A. In the short-wavelength limit A ^ the GHS 
therefore disappears leading to the standard ray dynamics of 
geometric optics. In the field of quantum chaos, the short- 
wavelength limit corresponds to the classical limit. For finite 
wavelength A the GHS can be regarded as a semiclassical cor- 
rection to the ray dynamics. 

The paper is organized as follows. The system is defined in 
Sec. Section Hn] presents the numerical analysis of local- 
ization of modes near ARCs. The semiclassical description in 
terms of the GHS is explained in Sec. |IV] Finally, a summary 
is given in SecIV] 



II. THE SYSTEM 

As model system we choose an optical microdisk with el- 
liptical cross section. Microdisks allow to confine light by 
total internal reflection at the boundary which is highly rel- 
evant for various research fields and applications 120]. Mi- 
crodisks with deformed cross section have attracted consider- 
able attention in the context of quantum chaos in open sys- 
tems I2II ll^l- Elliptical microdisks have been studied both 
theoretically ll23l l24ll and experimentally il25i l26ll . As semi- 
major axis we set a — R{1 + e) and as semiminor axis 
b = R/{1 + e). The area of the ellipse irab does not de- 
pend on the deformation parameter e. The eccentricity is 
given by e = ^1 — {b/aY = — 1/(1 + e:Y- We choose 
R = 1/im and the effective index of refraction n = 3.3 ac- 
cording to recent experiments on circular microdisks ifTTIl . 
Maxwell's equations for the transverse magnetic (TM) po- 
larized modes Ez{x,y,t) = ■ip{x,y)e~^'^'^ reduce to a two- 
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dimensional scalar wave equation |12 



(1) 



with frequency uj and the speed of light in vacuum c. The 
wave function tp and its normal derivative are continuous 
across the boundary of the cavity. At infinity, outgoing wave 
conditions are imposed. With these open boundary conditions 
the modes decay as light can leak out of the cavity. This re- 
sults in quasi-bound states with complex frequencies uj in the 
lower half-plane. The real part is the usual frequency. The 
imaginary part is related to the lifetime r — — l/[21m(a;)] 
and to the quality factor Q = — Re(a')/[2Im(ti;)]. 



III. AVOIDED RESONANCE CROSSINGS 

We compute the modes numerically using the boundary el- 
ement method |28]. Figure[T]shows the calculated normalized 
frequency fl = ojR/c = kR as function of the eccentricity e. 
The considered region around Re(ri) — 6.5 corresponds to 
a vacuum wavelength of about 970 nm. We observe that the 
real part of the frequency performs an avoided crossing, and 
the imaginary part shows a formation of a short-lived (large 
Im(f7)|) and a long-lived state (small |Im(fi)|). This is an 
ARC of external coupling in the strong coupling regime ifTsIl . 
Observing an avoided crossing is the more remarkable when 
recalling that the corresponding closed system with vanish- 
ing wave intensity on the boundary - the so-called elliptical 
billiard ||29[ [30ll - belongs to the class of integrable systems 
which typically do not show avoided level crossings [31]. 



6.55 



6.5. 



6.45 

-0.01 
-0.02 

% -0.03; 

-0.04 



iA"~-~^^ 









^ ^ 




:B 
















1 1 r.-< 







1 1 1 1 1 1 1 


1 1 1 1 i_ 










B ~ - - 










1 1 1 1 1 rrT*>Tri 1 1 







0.83 


0.84 0.85 
eccentricity 


0.86 


0.87 



FIG. 1 : An avoided resonance crossing in the elliptical microcavity. 
Plotted are the complex frequencies Q. as function of the eccentric- 
ity. The real part shows an avoided crossing, whereas the imaginary 
part shows the formation of a long-lived mode C and a short-lived 
mode D. Dots mark the frequency of the modes shown in Fig-O 

The corresponding spatial mode structures are depicted in 
Fig.|2] The modes A and B at the left hand side of the ARC 
and the modes E and F at its right hand side look very similar 
to the states in the elliptical billiard Ii29il . These states can be 



labelled by two mode numbers n\ and 712. With n\ (n^) we 
count the nodal lines in radial (azimuthal) direction in an ellip- 
tical coordinate system ll29ll . Note that, strictly speaking, the 
wave function of a quasi-bound state does not possess nodal 
lines, only the real and the imaginary part of the wave function 
does. For modes A and F we find (ni, = (2, 26), and for 
modes B and E we find {n\,ni) = (3, 18). 
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FIG. 2: (Color online) Calculated near field intensity of modes with 
the same labels as in Fig.[T] The long-lived "bowtie" mode C and the 
short-lived mode D show localization along periodic ray trajectories 
(dashed lines). 

The hybridized modes C and D, however, have no coun- 
terparts in the elliptical billiard. They are superpositions 
of the original modes A and B (or E and F) shown in 
Fig. |2] We can observe a clear accumulation of intensity 
on periodic ray trajectories. As these ray trajectories are 
marginally stable [ 29], such kind of modes have been referred 
to as "scarlike" in Ref. L15J. True scars refer to the local- 
ization of modes along unstable periodic ray trajectories in 
chaotic systems fvf]. It is customary to use the terminology 
"scars" also for systems with mixed phase space, see, e.g., 
Refs. lf32[ l33i l34[ l35[ [36[[ . Scarring has been observed also in 
microcavities "IMBqI liol |4ll| . 

What causes the lifetime enhancement of the bowtie-shaped 
mode C and the lifetime reduction of mode D7 In Ref. lUSll 
it was shown that the lifetime enhancement (reduction) of 
modes in rectangular cavities is due to destructive (construc- 
tive) interference at dielectric corners. Scattering at dielec- 
tric corners is the main decay channel in this kind of cavi- 
ties. For the elliptical cavity there are no corners. To under- 
stand the lifetime modifications in this and other cavities with- 
out corners it is convenient to analyse the (emerging) Husimi 
function ll42ll . representing the wave analogue of the phase 
space. A phase space representation of the ray dynamics is 
the so-called Poincare surface of section (SOS), as shown in 
Fig. [3] Whenever the trajectory leaves the cavity's bound- 
ary after being reflected, its position s (arclength coordinate 
along the circumference) and tangential momentum p — sin x 
(the angle of incidence x is measured from the surface nor- 
mal) is recorded. It is a special property of integrable systems 
that each trajectory is confined to a so-called invariant curve 
(a torus in the full phase space). The invariant curve is ex- 
emplarily indicated by the thin dots in Fig. [3] The symme- 
try properties of these curves in the SOS are determined by 
the discrete spatial symmetries of the billiard under investiga- 
tion and the p — > ~p symmetry. The latter is due to time- 
reversal symmetry and hard-wall boundary conditions. The 
region between the two critical lines for total internal reflec- 
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tion sinxc = ±l/n is called the leaky region. A ray which 
enters this region escapes according to Snell's and Fresnel's 
laws. 
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FIG. 3: (Color online) Emerging Husimi function (shaded regions) 
of the modes in Fig. (2] The horizontal lines are the critical lines 
sinXc ~ ±l/n. Thin dots mark some invariant curves in the phase 
space of the ray dynamics. Thick dots represent the periodic ray 
trajectories on which the modes C and D are localized. 

From Fig. 13 we can learn that far away from the ARC, the 
modes A, B, E, and F follow roughly the ray dynamics as 
they are located near invariant curves. Note that the p —p 
symmetry present in the corresponding billiard is broken in 
the open microcavity ll43l HJ]. By comparing modes A and 
B with modes C and D, it is apparent that near the ARC a 
rearrangement of the phase space structure takes place. The 
mode D has an enhanced intensity inside the leaky region 
leading to a reduction of its lifetime. The mode C, however, 
has strongly suppressed intensity in the leaky region of phase 
space which increases its lifetime. Moreover, the phase space 
picture shows clearly the localization of the modes along pe- 
riodic ray trajectories (thick dots in Fig.|3]l. 

In the following we demonstrate that the appearance of lo- 
calized modes at ARCs is not necessarily associated with the 
formation of short- and long-lived modes. Figure |4] shows 
an example. Instead of the formation of a long- and a short- 
lived mode we see a crossing in the imaginary parts which is 
a signature of an ARC of internal type in the strong coupling 
regime lfl4ll . The corresponding mode patterns are depicted in 
Fig. |5] The unperturbed modes A and B (E and F) show a 
whispering-gallery structure similar to modes in a circular mi- 
crodisk. The mode numbers are (rii, 712) = (1, 38) for mode 
A and (0, 46) for mode B. Also this type of ARC leads to 
hybridized modes, C and D, having strong localization along 
periodic ray trajectories. 

The localization is also clearly visible in the phase space 
representation in Fig. |6] Moreover, this representation ex- 
plains why there is no formation of short- and long-lived states 
in this case. Mode A lives closer to the leaky region than 




FIG. 4: An avoided resonance crossing with an avoided crossing in 
the real part of the frequency and a crossing in the imaginary part. 
Dots mark the frequency of the modes shown in Fig. [5] Inset shows 
a magnification of the avoided crossing. 




FIG. 5: (Color online) Calculated near field intensity of the modes 
with the same labels as in Fig. |4l The rectangular-shaped mode C 
and the diamond-shaped mode D show localization along periodic 
ray trajectories (dashed lines). 



mode B. Hence, the tail of mode A has a larger contribution 
inside the leaky region (even though the intensity is so weak 
that it cannot be seen in the figure). This explains why the life- 
time of mode A is much shorter (about a factor of 10) than the 
lifetime of mode B. Due to their very different contributions 
in the leaky region a superposition of modes A and B cannot 
lead to a significant cancellation of phase space intensity in 
the leaky region. Consequently, short- and long-lived modes 
cannot form. 

The observed localization patterns can be explained for the 
cavities where the corresponding billiard is integrable. The 
hybridised modes, e.g. C and D in Fig.|5] are superpositions 
of the unperturbed modes far away from the ARC, e.g. A and 
B in Fig. |5] As the eigenstates of an integrable billiard the 
unperturbed modes have a simple structure of nodal lines. 
Therefore they can be labelled by a pair of mode numbers 
{ni,n2) and (mi, 7712). The morphology of superpositions 
<^'4'ni,n2 + /3V'mi,m2 is Strongly influenced by the differences 
\ni — TOil and \n2 — m2|. If the differences are small num- 
bers compared to and rrii, then the mode pattern can be de- 
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FIG. 6: (Color online) Emerging Husimi function (shaded regions) of 
modes in Fig. [5] The horizontal lines mark the critical lines sin Xc = 
±l/n. Thin dots correspond to invariant curves in the phase space 
of the ray dynamics. Thick dots represent the periodic rays on which 
the modes C and D are localized. 



composed in a strongly oscillating part and a weakly varying 
envelope. The weakly varying envelope defines the localiza- 
tion pattern. Figure [T] illustrates this scenario for the simpler 
case of a rectangular billiard. The modes can be computed 
analytically: 

'^n^,ny{x,y) = siii(J^x^ sm(^^^y^ (2) 

if < a; < i? and < y < eR; otherwise i!n^,ny{x, y) = 0. 
The positive integers Ux, Uy count the number of nodal lines 
in X- and y-direction (boundaries a; = and y — are not 
counted), respectively. Figures |7|a) and (b) show the cases 



(10,7) and {mx,my) = (12,5). It can be seen 



that the superposition '0io,7~V'i2.5 has reduced intensity at the 
corner regions due to destructive interference. Correspond- 
ingly, the intensity near the centre points is enhanced due to 
constructive interference. The opposite is true for the superpo- 
sition tpioj + '!/'i2.5- In both cases the localization resembles 
scarring on periodic ray trajectories of short period. 



IV. GOOS-HANCHEN SHIFT 

In the previous section we demonstrated the existence of 
localized modes in optical microcavities near ARCs. But why 
is the light intensity in such a situation localized along peri- 
odic ray trajectories? To find an answer to this question we 
compare in Fig.[8]the incident and the emerging Husimi func- 
tion ll42ll forp = 0.75. The incident (emerging) Husimi func- 
tion represents the optical mode before (after) the reflection 
at the boundary of the microcavity. A lateral shift of about 
As = 0.011 •Smax can be observed. We interpret this shift as 
the Goos-Hanchen shift (GHS). 



4 ^ 




FIG. 7: (Color online) Eigenmodes and superpositions of eigen- 
modes in a rectangular billiard, (a) 11)10,7 {x, y), (b) '(/'i2,5(a;, y), (c) 
IV'10,7 — ■!/'i2,5p, and (d) IV'10,7 + V'la.sf- Lines mark periodic ray 
trajectories. 




FIG. 8: Incident (solid line) and emerging (dashed) Husimi function 
of mode C in Fig.|6]for p — sin x = 0.75. 



In the present section we develop a semiclassical descrip- 
tion based on the GHS. We restrict ourselves to systems whose 
closed counterpart is integrable. Examples are the ellipti- 
cal and the rectangular cavity. As well known in the field 
of quantum chaos, integrable closed systems typically show 
level crossings rather than avoided level crossings. In the fol- 
lowing we consider the crossing point of energy levels where 
the levels form a pair of degenerate levels. To do so, we use 
the quantum mechanical notation for the solutions of the wave 
equation. We demonstrate that the degeneracy of energy lev- 
els implies degenerate ray motion. Degenerate ray motion 
takes place on resonant tori which are foliated by periodic ray 
trajectories, in contrast to generic tori which support quasi- 
periodic motion. 

Let us consider two states each labelled by a pair of quan- 
tum numbers (ni, 722) and (mi, 7712), respectively. In the limit 
of large quantum numbers the semiclassical EBK quantization 
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rule m 



(a) 



S„,,„, = H[{m + ^)h,{n2 + ^)h], (3) 



Er, 



= H[{m, + ^)h,{m2 + ^)h] (4) 



/32 



is expected to give accurate results; H is the classical Hamil- 
tonian, and /3i are Maslov indices ll46ll . Note that for sys- 
tems with separatrices, Maslov indices can be defined only 
locally in phase space ll 29i WA l48ll . However, if both states 
belong to the same type of classical motion then a,; and Pi are 
uniquely defined. 

Moreover, we assume in the following that \ni — mi\ <C 
riijiTii. In this regime we can use the Taylor expansion in 
An," — nij — n," 



dH ^ ^ dH ^ ^ 
— — Anifi + — — An2ft . 
oh ah 



In the case of degenerate levels -Em^^ma — ^ni.n 
write Eq. (|5]i with the frequencies uoi — dH/dh as 

oji Aril + ti)2Ari2 — . 



(5) 



(6) 



Since Arii are integer numbers, Eq. (|6]l requires that the fre- 
quency ratio uji/u!2 is a rational number As a consequence 
the corresponding torus is resonant and composed of periodic 
ray trajectories. This line of arguments shows that the degen- 
eracy of quantum levels is related to the appearance of peri- 
odic ray trajectories for the case of integrable closed systems. 

When an integrable system is subject to a smooth and 
generic perturbation then on the quantum (wave) side the de- 
generacy of the levels is lifted. As a consequence, a level 
crossing turns into an avoided level crossing. In the case of 
microcavities the perturbation is given by the open boundary 
which allows for leakage of light out of the cavity. On the 
classical (ray) side, a generic perturbation breaks up the reso- 
nant tori into a chain of stable periodic orbits (ray trajectories) 
surrounded by islands and unstable periodic orbits surrounded 
by a chaotic layer This scenario is proven by the Poincare- 
Birkhoff theorem [49^ 50]. The fate of the nonresonant tori is 
determined by the Kolmogorov-Arnol'd-Moser (KAM) theo- 
rem lf5lll52ll53ll . It states basically that most of the nonreso- 
nant tori are conserved in the perturbed system. 

The basic idea here is to consider the GHS as a wavelength- 
dependent perturbation of the ray dynamics. The GHS is a 
direct consequence of the open boundary. The general sce- 
nario is as follows. The GHS creates stable and unstable pe- 
riodic orbits according to the Poincare-Birkhoff theorem. The 
size of the corresponding islands surrounding the stable orbits 
in phase space depends on the ratio of wavelength to cavity 
size. We get smaller islands for larger normalized frequen- 
cies n — kR and larger islands for smaller fi. Following 
the semiclassical eigenfunction hypothesis f54','55'] the islands 
can support modes. In addition, modes can be scarred along 
the unstable periodic ray trajectories. This semiclassical de- 
scription is summarized in Fig. |9] together with the wave de- 
scription. It is important to note that the nonintegrability of 
the augmented ray dynamics explains the existence of ARCs. 
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FIG. 9: (Color online) From the open boundary to localized modes 
near ARCs in systems with integrable closed counterpart, (a) Wave 
description and (b) semiclassical description with augmented ray dy- 



The GHS As for a TM polarized plane wave reflected at 
a planar dielectric interface is according to K. Artmann lf56ll 
given by 



kAs = Re 



smx 



\/ v? sir? X ' 



sin^ X 



(7) 



The angle of reflection is equal to the angle of incidence x- 
Figure [TO] shows that the r.h.s. of Eq. (jTjl has singularities at 
the critical angle sinxc = and at sinx 1- These sin- 
gularities disappear in the case of incoming waves with finite 
spatial extension. For Gaussian beams of width a Lai et al. 
derived an analytical expression [fs?']. However, their result 
is restricted to ka ^ 1, and exhibits artificial singularities 
in the low-fccr regime considered here. Therefore, we use the 
Artmann result dTji for long-lived ray trajectories which stay 
well within the interval sin x G [0.35, 0.9]. Note that the GHS 
As disappears in the semiclassical (short-wavelength) limit 
k ^ oo. The GHS at curved interfaces has been computed nu- 
merically in Refs. ll58il59|] where, away from the critical angle 
in the regime of total internal reflection, reasonable agreement 
with the Artmann result has been found. No analytical formu- 
las are available for the general case. 

To compute the augmented ray dynamics we combine 
the billiard map (sj+ijPj+i) — M{sj,pj) with the Goos- 
Hanchen map G(A) defined by 



Sj+i = Sj + As{pj, A) , 
Pj+i = Pj ■ 



(8) 
(9) 



As the map G(A) preserves the phase space area, the com- 
bined map MG{X) is Hamiltonian. In the short-wavelength 
limit we recover the standard ray dynamics, i.e. M ~ 
limA^o MG(A). 

Figure[TT]shows the phase space structure of the augmented 
ray dynamics for the situation of the rectangular mode at 
e K, 0.649. The left inset shows that the GHS creates small 
islands in phase space. The center of these islands correspond 
to stable periodic orbits. The phase space in Fig.[TT|has been 
computed with Eq. ^ for VL = kR = 8.2. It is important to 
mention that the overall island structure does not depend on 
fi; only the size of the islands does. The same island struc- 
ture is also obtained when applying the formula of Lai et al. 



FIG. 10: Goos-Hanchen shift As at a planar interface according to 
Eq. Q. The function As is antisymmetric with respect to x ~* —X' 
i.e. the shift measured along the direction of momentum is always 
non-negative. If|sinx| < 1/n then As = 0. The dashed line marks 
the critical line sin Xc = 1/n. 
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for Gaussian beams 115711 in the regime of large ka where the 
formula is valid. From Fig.[TT|it can be clearly seen that the 
Husimi distribution is correlated with the island chain. Note 
that for both the Husimi distribution and the ray dynamics 
the outgoing component is shown, i.e. we plot the emerg- 
ing Husimi distribution ll42ll and the orbit just after reflection. 
Figure [TT] reveals that the rectangular mode is not a scar with 
respect to the augmented ray dynamics but a mode living in 
a stable island chain following the semiclassical eigenfunc- 
tion hypothesis. While the center of the Husimi function and 
the corresponding island agree well, the orientation of the two 
objects is different. This can be understand by the way the 
Husimi distribution is constructed. The Husimi distribution 
can be defined as Wigner distribution (restricted to the bound- 
ary of the cavity) folded with Gaussian wave packets of min- 
imal uncertainty (coherent states) |42]. The relative uncer- 
tainty and the orientation of the wave packets can be chosen 
freely, but usually the orientation is chosen in accordance with 
the s- andp-axes of phase space, as it is done in Fig.[TTl How- 
ever, one can also choose the orientation in accordance with 
the orientation of the island under consideration. Such "tilted 
coherent states" have been used, e.g., in |60]. Moreover, note 
that the Husimi function for billiards and cavities itself is a 
semiclassical approximation |42]. 

The discussed ambiguity of the Husimi distribution is not 
present for the corresponding wave function in real space 
shown in the right inset of Fig. [TT] In this real space repre- 
sentation the GHS is directly visible. The real space plot of 
the mode and the stable periodic ray trajectory shows much 
better agreement than in the case of the standard ray dynam- 
ics; cf. Fig. |5] In fact the agreement is surprisingly good 
having in mind that we use the GHS for a planar interface and 
that Vl = 8.2 is deep in the wave regime. 

The strong leaking of the Husimi function out of the islands 
in Fig.[TT]can be confirmed by computing semiclassically the 
number of modes m fitting into an island of size A. From 
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FIG. 11: (Color online) Emerging Husimi function of rectangular 
mode (mode C in Figs. [5] and [6} compared to the Poincare SOS of 
the augmented ray dynamics above the critical line sinxc = 1/n. 
The thick dots mark the stable periodic ray trajectory at the centre of 
the islands (thin dots). Left inset: magnification of one of the islands. 
Right inset: mode structure and stable periodic ray trajectory in real 
space. Parameters are as in Figs.|5]and|6] 

EBK quantization for the action / = ^ <f pds = nih with 
p — nhk sin x one gets 



nfls„ 



2nR 



'-A{n) 



(10) 



where A{n) is the dimensionless phase space area spanned 
by sinx and s/smax; the frequency = kR is considered as 
real-valued for notational convenience. Numerically we find 
A ^ 1.53 X 10"'^ and m = 0.04. The value of m appears to 
be too small, however, it is well-known that also small islands 
can accommodate modes. This is usually explained by the 
confinement due to cantori in the phase space region around 
the islands |32. 61]. In our case the confinement is provided 
by the KAM tori (the unbroken tori of the elliptical billiard) 
surrounding the islands visible in the left inset of Fig.fTTI 

Figure [12] shows the same phase space as before but to- 
gether with the Husimi function of the diamond mode. In- 
terestingly, this mode is localized on an unstable periodic ray 
trajectory of the augmented ray dynamics. This periodic ray 
trajectory is suiTounded by a tiny chaotic layer which is not 
visible in Fig. [12] Hence, the diamond mode is a true scar 
in the augmented ray dynamics. Again, the real space plot in 
the inset shows much better agreement than in the case of the 
standard ray dynamics; cf. Fig. [5] 

Figure[T3]eventually depicts the case of the bowtie mode C 
with e « 0.845, cf. Figs. [2] and [3] Again, the GHS creates 
islands in phase space. The inset of Fig. [13] shows that the 
augmented ray dynamics fits slightly better the mode struc- 
ture than the standard ray dynamics shown in Fig. [2] In the 
case of Fig. [T3] we find for the area A — 4.7 x 10"'^ and 
the number of modes in the islands m — 0.11. The scenario 
for the hybridization partner, mode D, results in an unstable 
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FIG. 12: (Color online) Emerging Husimi function of diamond mode 
(mode D in Figs. |5] and ^ compared to the Poincare SOS of the 
augmented ray dynamics. The thick dots correspond the unstable 
periodic ray trajectory. Inset: mode structure and unstable periodic 
ray trajectory in real space. Parameters are as in Figs.|5]and|6] 



periodic ray trajectory (not shown). Here, the GHS affects 
only the overcritical reflections in the elongated region of the 
ellipse, whereas the shift is zero at the remaining reflection 
points because of normal incidence. 




0.4 0.6 



FIG. 13: (Color online) Emerging Husimi function of bowtie mode 
(mode C in Figs.|2]and[3]( compared to the Poincare SOS of the aug- 
mented ray dynamics. The thick dots marks the stable periodic ray 
trajectory at the centre of the islands. Inset shows the mode structure 
and the stable periodic ray trajectory in real space. Parameters are as 
in Figs.[2]and[3] 

Figure [T4f a) shows how the size of the islands scales with 
the frequency fi in the case of e « 0.845. For large enough 
fi we see a power law dependence with exponent k, —0.3. 
From Eq. ( fTOl i it then follows that the number of GHS modes 



increase as Jl"-^ in the semiclassical limit. Note that the first 
data point with il = 8.2 is not used in the fitting procedure 
since the corresponding island is strongly influenced by the 
critical line; see Fig. [13] For the case e « 0.649, Fig. [T4jb) 
shows a power law with exponent « —0.72. Here the number 
of modes in the semiclassical limit scales as f2"-^^. Interest- 
ingly, the exponent depends on eccentricity. In both cases the 
number of GHS-induced modes increases as the semiclassical 
limit is approached (or the wavelength is fixed and the cavity 
size is increased). 




FIG. 14: Size of the island A created by the GHS as function of Q. 
in log-log plot; (a) eccentricity e ~ 0.845. The line is a linear fit 
with slope ~ —0.3. (b) e ~ 0.649. The line is a linear fit with slope 
^ -0.72. 

GHS-induced modes have been already predicted for the 
case of microdome cavities ll62ll . However, the relation to 
ARCs has not been revealed in Ref . 16211 . 



V. SUMMARY 

We have analysed the phase space structure of localized 
modes formed near avoided resonance crossings. Using the 
Husimi function we are able to demonstrate the mechanism 
for the reduction (enhancement) of lifetimes due to construc- 
tive (destructive) interference in real space leading to a reduc- 
tion (enhancement) of intensity in the leaky region of phase 
space. We have developed a semiclassical extension of the ray 
model to describe the wave phenomenon of localized modes 
for systems with integrable closed counterpart. This aug- 
mented ray dynamics includes the Goos-Hanchen shift. As 
the Goos-Hanchen shift breaks the integrability it naturally 
explains both the appearance of avoided resonance crossings 
and the localization of modes. According to the Poincare- 
Birkhoff theorem the perturbation creates stable and unstable 
periodic ray trajectories along which the modes localize. This 
scenario shows that half of the "scarlike" modes of Ref. [Tsl 
can be regarded as true scars in the augmented ray dynamics. 
As example we discussed an optical microcavity with an el- 
liptical cross section, but we expect that our results apply to 
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other cavity geometries as well. Moreover, we believe that the 
concept of an augmented ray dynamics can also be a ppli ed to 
quantum systems, in particular to soft-wall billiards ll32ll63ll . 
in order to semiclassically describe quantum effects in a con- 
ceptually simple way. 
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